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Abstract 

We consider the so called Moran process with frequency dependent fitness 
given by a certain pay-off matrix. For finite populations, we show that the final 
state must be homogeneous, and show how to compute the fixation probabili- 
ties. Next, we consider the infinite population limit, and discuss the appropriate 
scalings for the drift-diffusion limit. In this degenerated parabolic PDE 

is formally obtained that, in the special case of frequency independent fitness, 
recovers the celebrated Kimura equation in population genetics. We then show 
that the corresponding initial value problem is well posed and that the discrete 
model converges to the PDE model as the population size goes to infinity. We 
also study some game-theoretic aspects of the dynamics and characterize the 
best strategies, in an appropriate sense. 

1 Introduction 

Since the beginning of modern evolutionary theory, the study of the dynamics of a 
mutant gene in a population has attracted attention |TQ1 HU E3 EZ1 EE] ■ It has been 
known for a long time that a mutant gene will be, eventually, either fixed or lost. The 
final result depends not only on natural selection but also on chance j2Uj . 

The most natural attempt to describe mathematically the evolution of a mutant 
gene uses a discrete model for a finite population. The question of finding a consistent 
model for the infinite population is then a natural one. This is called in the physical 
literature the "thermodynamical limit", and it is a classical subject on that field. See, 
for example, [H|. 
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When we consider the infinite population limit, it is natural to have continuous 
variables where we previously had discrete ones. For example, if N denotes the size of 
the population, the possible fractions of mutants are 0, 1/N, 2/N, In the limit 

N — > oo this fills the entire real interval [0, 1]. Time should be rescaled accordingly, 
such that that the probability of fixation, for a given fraction x over a given time span 
t, does not depend significantly on the size of the population. In the infinite limit, we 
obtain a partial differential equation (PDE). This PDE is an approximation for large 
iV of the discrete process, and as such must present diffusion to the boundaries as a 
continuous representation of the fact that the mutant gene will be eventually fixed or 
lost. 

A different approach to the same problem is to consider continuous models from 
the beginning. This leads to two distinct modeling paradigms: the first one uses 
ordinary differential equations (ODEs) to model the evolution of the fraction of mu- 
tants. The most widely used equation in this context is the replicator dynamics and 
its variations ^B] • The second one, which will be further developed in this work, uses 
PDEs to model the evolutionary process. This approach goes back, at least, to the 
seminal work by Kimura |20j . where it was used to model the diffusion of mutant 
genes. More explicitly, Kimura considered the probability of fixation of a mutant 
gene that in a given time is present in a certain fraction of the population. Here, we 
will deduce the Kimura equation as a particular case of our work, where a PDE will 
be obtained from the more basic discrete process, in the infinite population limit, for 
the diffusion of mutant genes. 

In this work, we consider a simple evolutionary process, called the Moran process, 
introduced in j2H] (used, e.g., for cancer dynamics [T3 I2H], paleontology [23], phy- 
logeny j21j, genealogy [0], and epidemiology 36J) for a finite population and obtain 
a partial differential equation as its thermodynamical limit. 

Our starting point is evolutionary game theory. We consider a finite population of 
fully connected interacting individuals through a certain pay-off matrix. We start by 
proving that for any finite population (of size N) of two-types, one of the types will 
be fixed after long enough time. The thermodynamical limit is then obtained as a 
PDE that approximates the finite population dynamics for large N. We consider two 
different scalings for the time-step At, namely: At = 1/N (drift limit) and At = 1/N 2 
(drift-diffusion or simply diffusion limit). In the second case it is also important to 
introduce the so called weak selection limit (pay-offs go to 1, when population goes 
to infinity). We also show that the most interesting equations appear in the drift- 
diffusion limit. 

All the equations found in the limit are degenerate, i.e., the diffusion coefficient 
vanishes on the boundaries. The mathematical theory for such equations is not as 
well developed as for the non-degenerate case. There are the classical books [21 Ej- In 
particular, [2| proves existence and uniqueness for the equation obtained by Kimura. 
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For more recent works, see also jUEj. 

If we impose no diffusion in the PDE model, the solution can be decomposed in 
point dynamics, where each fraction evolves through the replicator dynamics. The 
stationary states and long time behavior of the replicator dynamics are, however, 
different to the ones obtained as the thermodynamical limit of the final states of 
discrete populations, showing that the diffusion is essential to understand the discrete 
dynamics. 

The PDE model allows the introduction of a relation of dominance between two 
different strategists that turns out to be, in its dynamical features, identical to the 
flow of the replicator dynamics. We also show that the best possible strategy in the 
finite, but large, population case is given by the evolutionarily stable strategy (ESS) 
of the game [THl I32j. This clarifies the relation between a homogeneous population 
playing mixed strategies with given frequencies and a mixed population, with constant 
fractions, playing pure strategies, i.e., the difference between evolutionarily stable 
strategies and evolutionarily stable states. 

If the fitness for the individuals in the population is frequency independent, the 
resulting equation is equivalent to a well-know equation of population genetics, in- 
troduced by Kimura [20], describing the probability of fixation of a mutant (with 
frequency independent fitness) in a population. It turns out, that the equation de- 
rived in this work and the one introduced by Kimura are a forward/backward pair of 
equations. 

It is important to note that the perception that the Moran process (at least in 
the frequency independent case) is related to diffusion process is not new jU E] . The 
compatibility between finite populations simulations and the ESS, defined in the 
continuous case, are also studied in [T21 EH I2H] • 

The structure of this work is the following: In Section |21 we introduce the (finite 
population) Moran process and study its properties. In particular, we prove that the 
final state will be always homogeneous. In Section El we introduce the drift-diffusion 
scaling and obtain a PDE as the thermodynamical limit of the Moran process. We also 
study its dynamic features from the strategic point of view. In Section 01 we consider 
the no diffusion case and compare the PDE obtained with the replicator dynamics. 
In Section we particularize all results to the frequency independent case and in 
Section El we study the drift scaling. Finally, in Section [7| we point new directions for 
this work, showing how the tools developed here can be applied to different dynamics. 

2 The frequency dependent discrete case 

We consider a fixed size population with two types of individuals: A and B, say. At 
fixed time steps, we choose one of the individuals to be eliminated at random and 
replace it by a newborn which can be of either type. This newborn is obtained as a 
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Figure 1: The Moran process: from a two-types population (a) we chose one at 
random to kill (b) and a second to copy an paste in the place left by the first, this 
time proportional to the fitness. 



copy of one of the remaining individuals with probability proportional to its fitness. 
See Figure Q for an illustration. This process is called the Moran process [23] • 

Let P(t, n, N) be the probability that there are n type A individuals at time t in a 
population of fixed size N. We define c+(n, N) (co(n,N) and c-(n,N), respectively) 
as the probability (independent of time) that the number of mutants changes in time 
t from n to n + 1 (to n and to n — 1 respectively) in time t + At. We assume that 
these transition probabilities are proportional to the fitness 4>a and 4>b of types A 
and B respectively; thus we have: 

c+(n,N) = —— 7T— , 1 

N ncpA + [N — n — l)<p B 

n (n-l)<j) A N-n (N - n - l)(f> B 

C ° ln ' ] N (n — \)<j>A + {N- n)4> B N n0 A + (N - n - \)<f> B ' 

f at\ n {N-n)<f>B /q\ 

C ~ {n ' N) = N { n-l)4, A + {N-n)4, B - (3) 

From that, we may easily write an equation for the evolution of P: 

P(t + At,n,N) = c + (n-l,N)P(t,n-l,N) + co(n,N)P(t,n,N) 

+c_(n+l,N)P(t,n + l,N) . (4) 

After imposing the boundary conditions P(t, —1,N) = P(t,N + 1,N) = 0, Vt > 0, 
we conclude that the previous recursion is valid for t > and n = 0, 1, • • • ,N. 

Originally, the Moran process was defined with a frequency-independent fitness, 
i.e., 4>a,b were independent of the particular composition of the population. We con- 
sider, however, the frequency-dependent case, and we obtain the results for frequency- 
independent populations special case. 

Now, we obtain the fitness. For that, we first consider a two players game, with 
pay-off matrix given by: 
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where I and II are two pure strategies and A,B,C,D > 0. We call an ^-strategist 
an individual that plays I with probability q and II with probability 1 — q. 

We assume that the two types play two (possibly) different strategies, E qi and 
E q2 . The pay-off matrix is then given by 



where 









E qi 


E Q2 










E qi 


A 


B , 












C 


D 






A 


:= q*A + qi (l- 


<?i)(£N 


-c) + 


(l-qi) 2 D, 






B 


:= q^A + q^l 


- q 2 )B + (1 - 


qi )q 2 C + (l 




- q 2 )D , 


C 


:= qi q 2 A + (l- 


qi)q 2 B + qi(l 


- q 2 )C+{l 




- q 2 )D , 


b 


:= qlA + q 2 (l- 


q 2 )(B^ 


-c) + 


(1 - q 2 ) 2 D . 







(5) 
(6) 
(7) 
(8) 



For simplicity, we consider in this section only pure strategists, i.e., E\- and Eq- 
strategists for type A and type IB individuals respectively. The general case follows 
easily from the results in this section replacing (A, B, C, D) by (A, B, C, D). 

We identify fitnesses and pay-offs, and then we have that the fitnesses for I- and 
Il-strategists, for a population with n I-strategists, are given by 



n—\. N — n „ 

A+— — -B , n = !,-•• ,JV 



N-l N-l 

n „ N — n — 1 



N — 1 . 



(9) 
(10) 



Then, the evolution iteration is given by Equation (@J) with transition coeffi- 
cients ©-© and P- finjl . 



2.1 The discrete dynamics 

A natural question is what are the steady states of the iteration defined by the Moran 
process. Here we show that the discrete model cannot have a non-pure equilibrium. 
Let us define the relative fitness as 



, , Mn) (A-B)n + BN-A ^ 
pN[n) = — = (C-D)n + (N-1)D > °' 
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Also, let 



n ( N-n \ N-\ + {p-l)n 
fN{n) = — - and g N (n,p) = 



N \ N J "V 

Then it is a straightforward computation to verify that 

fN(n)p N (n) f N (n 



c+(n,N) 



g N (n,p N {n)) 



c_(n,N) 



g N (n - l,p N (n)) 



co(n,N) = l-f N (n){ PN ^', 4 



and 



g N (n,p N (n)) g N {n - 1, p N (n))J ' 

Let M be the iteration matrix of (J1J. Then MisaiV+lxiV+1, tridiagonal matrix, 
with entries given by 



Mii = Co{i,N), i = 0,...,N, 



(i, iV) and M i(m) = c_(z + 1, N), i = 0, . . . ,N - 1. 



From this, and the fact that pAr(n) > 0, it is easy to see that M is a nonnegative 
matrix. Since co(n, N) + c+(n, N) + c-(n,N) = 1, M is column stochastic. 

The answer to question raised in the beginning of this section is given by the 
following result: 



Proposition 1. Let M be as above and let P(t) = (P(t, 0), P(t, 1) 
Then 



1. 



lim M fc 

k— +oo 



/l 





\0 Fi 







'JV-l 





V 



where the F n satisfy 



F n = c+(n, iV)F n+ i + c_(n, iV)F n _i + c (n, N)F n , 
F = and F N = 1. 



P(t,JV)) 1 
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2. If 1 denotes the vector (1,1,..., l)t ; F = (F , F±, . . . , FnY and if (•, •, ) denotes 
the usual inner product, then we have that 

(P(t),l) = (P(0),l> and (P(t),F) = (P(0),F). 

In particular, the I 1 -norm of a nonnegative initial condition is preserved. 

Proof. For part 1, see the proof at appendix 1X1 

As for part 2, we first observe that, if a vector V satisfies M^V = V, then we 
have that 

(P(t + At), V) = (MP(t),V) = (P(t),MtV> = (P(t),V). 

Hence 

(P(t),V) = (P(0),V). 
From the fact that M is column stochastic, we easily conclude that 

M f l = 1, 

and the first invariant follows. For the second invariant, we observe that Equation 
(jll|) can be written in matrix notation as 



M f F = F, 

which concludes the proof. □ 

Remark 1. The two invariants described in part 2 of the proposition^ are the only 
invariants of the Moran process and play an important role in the determination of 
the correct continuous solution. 

Thus, the equilibrium states must have their mass concentrated in the extremes. 
The F n turns out to be the fixation probability of I-strategists, when the process start 
with n I-strategists. 

From the definitions of c*(n, N), we see that F n satisfies: 



n, T (n\F , ( n*r(n\ I 9N(n,p N (n)) \ p g N (n,p N {n)) 

PN {n)* n+1 - ^p N {n) + gjv(n _ liPJv(n)) )t n + gN{n _ ltPN{n 
F = and F N = 1. 

Equation ()12j) can be solved by writing 



-F i = n 

)) n_1 ' (12) 



H(n)= ^(n, PN (n)) ^ = _ 

pN{n)g N {n - I, p N (n)) 



s 



Then, ignoring the boundary conditions for the moment, we have that 

G n +i = H(n)G n , 

with solution given by 

n-l 



G n = G 1 Y[H{i). 



8=1 



Since 



n-l 



F n -F n . 1 = G 1 l[H(i), 



i=l 



we obtain, after applying Fn = 1 and F = 0, that: 

n k—1 



F n = G 1 J2U H ^ 

k=l i=l 

tJV fe-1 \ _1 

fc=l i=l / 

The expression given by ()13|) does not appear to yield a simple formula in the general 
case. However, compare the formulas found in Section where we study the case 
when the relative fitness is constant with respect to n. 

Remark 2. The coefficients obtained in the above analysis are for a Death/Birth 
process. For a Birth/Death process, they are simpler and are given by 

c + {n,N) . f»WP»W 



c_(n, JV) 



g N {n,p N {n)) 

fN(n) 
g N (n,p N (n)) 



where 



co(n,N) = 1- _ / N{n \ (l + p N (n)) 
gN{n,p N {n)) 



9N{n,p) = - 



N 

Also, in this case H(n) simplifies to 

H(n) ' 



p N (n)' 



Figure 2: Fixation probabilities for N = 20, A = 2, 5 = 1, C = 3 and D = 1. 
The points are taken from M 10000 , while the lines are obtained by numerically solving 



2.2 Numerical Results 

We numerically computed the M 10000 , for N = 20 and various relative fitnesses. The 
entries predicted to be zero by Proposition [T] were found to have magnitude less than 

io- 50 . 

Also, from these calculations, we extracted the fixation probabilities and compared 
them with the ones obtained by evaluating (JT3J) numerically. The result for a specific 
choice of fitness is displayed in Figure El For the case of frequency independent fitness, 
we can obtain explicit formulas for the fixation probability — see Section |S] — and we 
also compare with the fixation probabilities extracted from M 10000 in Figure 01 



3 The t her mo dynamical limit 

The aim of this section is to derive a continuous approximation, i.e., a PDE model 
for the discrete process described in the previous section. 

We define the probability density that at time t we have a fraction x 6 [0, 1] of 
type A individuals 

Pit xN N) n 
V(t,x,N):= V ' 1 ' = NP(t, xN, N) , with x = — , n = 0, 1, 2, . . . , N. 

N 

Furthermore, we assume that in the limit iV — > oo, V(t, x, N) converges in some sense 
to a function p(t, x) which is sufficiently smooth so that 

p(t,x±j^j= p(t, x) ± ±d x p(t, x) + ^d 2 xP (t, x) + 0(N~ 3 ) (14) 
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Figure 3: Fixation probabilities for constant r = C/A = D/B, when N = 20 com- 
puted from M 10000 together with the analytical fixation plotted as continuous func- 
tions of n/N; (a) r = 1 (b) r = 1.5. 

and re-write equation (jlj) to second order in N^ 1 as 



pit + At, x) — pit, x) 



N 



c { l ] + eg) + c W )p-( c^> - c^> ) d x p 



,(0) JO) 



(15) 



1 



N 3 



{c^~c^)d x p + l -{cf + c^)dlp 



where c* = c*\x), * = +, 0, — , i = 0, 1, 2, are defined by 



c (xN,N) 



c_( [ X +L)N,N 



C -(--i,iv) = cL 0) + icL 1) + ^. 



(16) 
(17) 
(18) 



Then 



+ 4 1} + c w = (Ax 2 + D(l 



ci 0) - cL 0) 



xf + (B + C)x{\ - x)Y ■ [A(A - C)x 4 
+ [B(B - D) + C(C - A) + 2C(5 - £>)) ,x 2 (l - x) 2 
+£>(£> -B)(l- x) 4 + 2ar(l - x) (A(B - D)x 2 -{A- C)D(1 - x))] 
x{\ - x) (x(A - C) + (1 - x)(B - D)) 



(19) 



Ax 2 + D(l - x) 2 + (B + C)x(l - x) 



(20) 
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If we impose that 



lim (A, B, C, D) = (1,1, 1, 1) , (21) 
lim N(A - 1, B - 1, C - 1, D - 1) = (a, b, c, d) , (22) 

JV— >oo 



we find 



lim (cf + 4 2) + c (2) ) = -4 



N 

I 1 1 1 l / / 



lim ( c2 } - c W ) = -2 + 4x 



N^oo 



and, from (liyfEZUl) . we have 
lim 

AT— >O0 



lim (c4 0) +cL 0) ) = 2x(l-x) 



lim iV (c+ ] + 4 1} + c a)N ) = -3a; 2 (a-6-c + rf)-2x(a-c-2(6-d)) + (d-6) , 
\im N (cf ] -c^) = x(l -x)(x(a-c) + (1 - x)(b - d)) . 

Finally, we divide Equation (JTSJl by At = N~ 2 (diffusive scaling), and take the limit 
N — > oo, to obtain 

d t p = [3x 2 (a — b — c + d) — 2x(a — c — 2(6 — d)) — (b — d)] p 
—x(l — x)(x(a — c) + (1 — x)(b — d))d x p 
+(-2)p + 2(1 - 2x)d x p + x(l - x)d 2 x p 

i.e., 

d t p = d 2 [x(l — x)p] — d x [x(l — x)(xa + (1 — x)f3)p] . (23) 

where a = a — c and (3 = b — d. We also define rj = a — (3. 

Supplementing Equation (J23j) we have the following conservation laws: 

d f 1 d f 1 

— / p(t, x)dx = and — / ip(x)p(t, x)dx = 0, 
dt Jo dt Jo 

where ^(x) is given in Theorem EJ 

Remark 3. It is important to stress that if we do not impose conditions 

there are another possible scalings. More precisely, if < PJ7j) still holds but / Tl%)) is 

replaced by 

lim N V (A — 1, B — 1, C — 1, D — I) — (a,b, c, d), < v < 1, 

t/ien another possible scaling is given by taking At = (l/N) 1+u and, in this case, 
we obtain \2ty) without the diffusion term. This equation is discussed in Section ^| 
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Moreover, if we drop \21\) - K2lfy . and only require that the payoffs have a finite limit 
when N goes to infinity, then yet another scaling is given by At = 1/N and, in this 
case, the equation for the probability density is given by 



d t p = d x 



x(l - x) (x(A - C) + (1 - x)(B - D)) 

p 



_x 2 (A - B -C + D) + x(B + C -2D) + D 
We analyze this equation in Section^ 



(24) 



Equation (j23j) is not readily covered by the usual theory of parabolic PDEs. How- 
ever, the analysis can be extended to obtain the following result: 

Theorem 1. 1. For a given p° G L 1 ([0,1]) ; there exists a unique solution p = 
p(t,x) to Equation P^| ) of class C°° (R + x (0, 1)) that satisfies p(0,x) = p°(x). 

2. The solution can be written as 

p(t, x) = q(t, x) + a(t)5 + b(t)8i, 

where q G x [0, 1]) satisfies without boundary conditions, and we 

also have 

a(t) = / g(s,0)ds and b(t) = / q(s, l)ds. 
Jo Jo 

In particular, we have that p G C ,00 (M + x (0, 1)). 

3. We also have that 

lim q(t,x) = (uniformly), lim a(t) = ir [p ] and lim b(t) = ni[p°], 

t— »oo t— »oo ' t^oo 

where tt and 7Ti are computed in Theorem [B Note that this means that the 
solution solution will 'die out' in the interior and only the Dirac masses in the 
extremities will survive. 

4- Assume p° G L 2 ([0, 1]) and let J{t) = J* x(l — x)q 2 (t, x)dx. Then, we have that 

J(t) < J(0)e~ 2Ao ', A > 0. 



See the proof at Appendix [Bj 

For completeness we show various numerical simulations for computing p(t,x). 
Due to display convenience we plot P(t,x) = (Ax)p(t,x), instead of p(t,x). See 
Figures iHIl 

We observe that p'(t,x) = p(t, 1 — x) also satisfies (J23|) changing the parameters 
(a,j3) — > (— (3, — a). Hence each computation actually yields solution for two set of 
parameters, just by reflecting the solution around the axis x = 1/2. 
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Figure 4: Solutions for P(t, x) for various times, when /3 — a — 0. This is the 
pure diffusive constant fitness case. Note the diffusion to the boundaries. The initial 
condition is given by p°(x) = Si/2(x). 




Figure 5: Solutions for P(0.1, x) for various values (3 and r\ := a — (3. Here, the initial 
condition is the same as in Figure 0] 



Figure 6: Solutions for P(t, x) when j3 = 2 and rj = for various times. This is the case 
of some drift with constant fitness. The initial condition is p (x) = 20a; 3 (1 — x), which 
is asymmetric with a peak at x = 3/4. Notice that the form of the initial condition 
together with the drift sign leads to a very rapid convergence to the equilibrium state. 
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Figure 7: Solutions for P(t,x) for various times, when (3 = 2 and 77 = 1. The initial 
condition is the same as in in figure El Notice that there is little difference from the 
computation with 77 = thanks to the form of the initial condition and to the order 
one size of the parameters. 



Figure 8: Same as Figure but with (3 
case. 



1 and 7] = 2. Same remarks apply in this 



Figure 9: Solutions for P(t, x) for various times when (3 = 10 and r\ = 20 with the 
same initial condition as in Figure |H1 The convergence for the equilibrium state is 
very fast also in this case. 
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Figure 10: Solutions for P(t,x) for various times when (3 = —20 and r] = —40. In 
this case the drift forces the solution to accumulate in the opposite direction of the 
initial large concentration. 




Figure 11: Solutions for P(t,x) for various times when (3 = 20 and r) = —40. Here, 
the convective term vanishes at x — 1/2. The effect is that, at first, the solution 
convected until the its peak reaches x — 1/2. Then it essentially stays there, while 
diffusion enforces the transport to the boundaries. In the second figure, the very ends 
of the interval are omitted for better view of the behavior in interior. 
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Figure 12: Solutions for P(t, x) for various times, when (3 = —20 and rj = 40, with 
initial condition p°(x) = 6x(l — x). In this case, the sign of rj drives the solution out 
of x = 1/2 to the extremes. The initial condition was chosen to be symmetric in this 
case to highlight this behavior. Also, as in the previous example, the second figure 
have the very ends of the interval omitted for a better view of the inner behavior. 

Notice that, for a — (3, (3 3> 1, we expect a behavior drift-dominated for interme- 
diate times. This means that, if x* — (3/{(3 — a) ^ (0, 1), then the solution will be 
convected until it reaches one of the boundaries, and then will diffuse to the steady 
state. Otherwise, depending on the sign of 7] the solution will either first concentrate 
near x*, and then diffuses to the boundary, or depart from x* in both directions 
towards the ends. Notice also, that the solutions are never smooth at the ends. Com- 
putations with different values of a and (3 produces qualitatively similar graphics. 

Now, let us go back to the general case, i.e., for E qi - and .^-strategists, instead 
of only pure strategists. Then, in a straightforward way (see also Equations ©-©), 
we define 

a := qfa + qi(l - qi)(b + c) + (1 - qifd , 
b := q 1 q 2 a + qt(l - q 2 )b + (1 - qi)q 2 c + (1 - gi)(l - q 2 )d , 
c ■= qi q 2 a + (1 - q 1 )q 2 b + q x {\ - q 2 )c+ (1 - gi)(l - q 2 )d , 
d := q la + q 2 {l-q 2 ){b + c) + {l-q 2 ) 2 d . 

Then, the equation for p, the fraction of E qi -strategists in the population is given 

I'.v 

d t p = d 2 x (x(l — x)p) — d x (x(l — x){x{d — c) + (1 — x)(b — d))p^j 

= d 2 x (x(l - x)p) - d x (x(l-x){x(a + {l-x)j3))p\ , (25) 

where a := a — c— (q 1 -q 2 )(qia + (l-qi)(3) and (3 := b-d = (qi-q 2 )(q 2 a + (l-q 2 )/3). 
Note that a — (3 = (qi — q 2 ) 2 (a — (3). Then, if q\ ^ q 2 and a ^ (3, then a ^ (3. 
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Theorem 2. For p(0, •) = p° G i| fl L°°([0, 1]), i/ie solution of Equation 
unique, non-negative, and accumulates on the boundaries, i.e., p°° := lim^oop = 
^ob ]^ + ^lb ]^; where vr [p ] = 1 — 7Ti[p ] and i/ie fixation probability of E qi 
strategists is given by 



£ Jy P°( x ) dx exp (-y 2 (qi - g 2 ) 2 ^/ - y(qi ~ <?2)(<?2a + (1 - 92)^)) dy 

Jo exp {-y 2 {qi - g 2 ) 2 ^ - - ^fea + (1 - g 2 )^)) dy 
Jo Jq^p ^) ex P (~i/ 2 ( q i - g2) 2 ^ - - g2)(g2« + (1 - 92)/?)) dy dx 
Jo ex P (-y 2 (<?i - <?2) 2 ^ - y(qi - 42) tea + (l - ? 2 )^)) dy 

Proof. It is enough to prove for gi = 1 and g 2 = (i.e., for Equation Q'23|)) and then 
change the result from (a,/3) to (&, (3). Existence, non-negativeness and convergence 
to the boundaries follows from Theorem ^ 

To obtain values 7Tj[p ], % = 1,2, we multiply Equation (J2l?j) by an d integrate 
from to 1. On assuming that p is such that integration by parts can be performed 
and that no boundary terms arise, we obtain that 



d t / p(t, x)ijj(x)dx 



x(l - x)p(t, x) O" 0) + ( x ( a ~ P)+ PW{x)) da; . 



Conservation laws are obtained solving i])"(x) + (x(a — (3) + (3)ip\x) = 0. Solutions 
are given by ip = cte (conservation of probability) and 



il>(x) 



rx i> i ~y 2 ~~Tr~ -vP ) dy, 



Remark 4. Notice that ip(x) is the continuous counterpart to the discrete fixation 
probabilities. 

Using that 



p°(x)ijj(x)dx = / (ir [p°]5o + iri\p°]5i) ijj(x)dx 



we get 



Jo IyP°( x ) dx ex P {-y 2 ^r - yf 3 ) d y 



J^exp (-y^-yp) d y 



fjJ^opl^-jU) dyd 



Jo ex P (-y 29 ^-yP) d y 

Finally, we change from a, f3 to a, f3. 



□ 
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Corollary 1. Ifp° = S x o, then 

/;° exp {-y\ qi - g 2 ) 2 ^ - yfa - q 2 )(q 2 a + (1 - q 2 )(3)) dy 



-Kl[5 x o] 



Jo exp (-y 2 (qi - g 2 ) 2 ^/ - y(qi - <?2)(g2« + (1 - 92)/?)) dy 



Definition 1. We say that E q2 dominates E qi (E q2 y E qi ) if, for any initial condition 
p° G L\ fl L°°([0, 1]), the probability of fixation for the strategy E qi is smaller than 
the one for the neutral case (the case qi = q 2 ), i.e., 



-l 

ni[p°} < ^[p°] := / xp°(x)dx . 
Jo 

We also say that E q2 S-dominates E qi if the above formula is valid for all p° = 6~ x o, 
x° e (0,1), i.e., 

tt^o] = J iy q2)i ^ y < x° Vx°G(0,l), (26) 

Jo F (n,<i2){y)dy 

where we defined the auxiliary function 

F (qi,q2)(y) ■= ex p(^-y 2 (Qi~Q2) 2 ^-^--y(qi-q2)(q2a + (i-q 2 )p)^j . (27) 

The following lemma shows that the two definitions above are in fact equivalent: 
Lemma 1. E q2 5-dominates E qi if and only if E q2 y E qi . 

Proof. We only need to prove the only if case. Let us consider any initial condition 
given by p° e L\ n L°°([0, 1]). Then 



Iq 1 Jy P°( x ) F (n,Q2)(y) dxd y 



Trib ] 



io J y 



Jo F {nm)Wy 
Jo fy Jjp°( z ) s ( z - x)F {q ^ q2) (y)dzdxdy 

Jo F (<n,<i2){y)dy 

— x)F( qim )(y)dydzdx 



Jo F (n,<i2)(y) d y 

] JlZisn 1 

Jo 1 F (m,c2)(y) d y 



Jo Jo L F (cn „ r! )(y)dy 



Now, we use Equation (J2S|) and conclude that 

7i"i[p ] < / / p°(z)8(z — x)xdzdx = / p°(x)xdx = 7rf [p°] . 



JO 



□ 



19 







r'/ n, if PlTin otiIaa if 


a > (3 > 


q* < 


g2 > gi 


a > > (3 


<f e (0, l) 


g2 < gi < g* or g 2 > gi > g* 


> a > (3 


g* > 1 


g2 > gi 


> (3 > a 


g* < 


g2 < gi 


(3 > > a 


<f e (o, i) 


gi < g2 < g* or q 1 > q 2 > q* 


(3 > a > 


g* > 1 


Q2 > gi- 



Table 1: Dominance relations for the non-degenerated (a ^ /3 ^ ^ a) thermody- 
namical limit of the frequency-independent Moran process, given by Equation (|25jl . 
with q* = 13/ \(3 -a). 



In view of this lemma, from now on, we consider only initial conditions of 5-type, 
i.e., p° = 5 x o. In order to prove dominance relations, we prove first the following: 

Lemma 2. // F< qim \ is increasing in the interval [0, 1], then E q2 y E qi . 

Proof. For p° = 5 x o, Equation (|2*U|) can be re- written as 

{y)dy < / (y)dy, Vx°G(0,l). 

Jo 

This equation can be interpreted as saying that the average of the function F^ qim ) in 
any interval [0, x° G (0, 1) is less than the average in the interval [0, 1], which is 
true whenever the function is increasing. □ 

Finally we prove the full relations of dominance for a 2 x 2 game. 

Theorem 3. Let E qi and E q2 , q±, q 2 G [0, 1], be two strategists in a 2 x 2 game, and 
let q* = (3/ {(3 — a). Then the relation of dominance is given by Tabled 

Proof. The proof consists in a long and tedious calculation proving that, for each 
range in Tabled the function F^ qim ^ is increasing. Then we use Lemma |5J □ 

The following corollary shows that the strategy E q * is the best possible strategy 
if (3 > > a. 

Corollary 2. If/3>0>a, then E q * y E q , Vg ^ q* := (3/{/3 - a) G (0, 1). 
Proof. For g 2 = g*, simplifies for 

F (g,q*)(y) = exp (~y 2 (q - q*) 2 ^-^j ■ 
For a — (3 < 0, this is an increasing function of y and this proves the corollary. □ 
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i 



Figure 13: Relation of dominance between E qi - and .^-strategists for given param- 
eters. Here, q* = (3/ {(3 — a). The first 2 figures (above) show dominance from pure 
strategies, the third one (a > > /3, below, left) shows that the pure strategies dom- 
inates their neighbors (and everybody dominates E q *) and the last one {(3 > > a, 
below, right) shows that E q * dominates any strategy. The arrow points from the 
dominated to the dominant. 

In order to finish the full picture of dominance, we need also the following: 
Lemma 3. If E q2 y E qi , then E qi )/- E q2 . 
Proof. First, we see, from Equation (|27j) . that 

F (qi,q 2 )(y) = F (q 2 ,qi) i 1 ~ V) F {qim) W • 

Then, we write 



1 f x 1 

- J o F (^)Wy = - G 



1 l-l-X 

F {q2,qi) 

{y)dy- / 
o Jo 



-^(gi,'?2)( 1 ) ■ 



Furthermore, 



F (qi,q2) {y)dy = / F {q2,qi){y)&v F {qim)( l ) ■ 

Jo 



Using the fact that E q2 y E qi , we find 



x 







F {Q2,qi){y)dy 



o 



pl-x 




/ F (q2,qi)(y) 




/o 





We re-arrange the terms and conclude that 



, „ / F (q2,qi) (y)dy > / F (q2:qi) (y)dy 
L x Jo Jo 



□ 
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Table Q together with Lemma El can be summarized in Figure IT31 
It is important to note also that in some references (see, e.g., [33]) it is said that 
selection favors strategy II replacing strategy I (in this case, we say that strategy II 
weakly dominates strategy I), in a finite population of size N, if a single type II mutant 
has fixation probability larger than 1/N, the neutral probability. Unfortunately, no 
sound generalization of this concept can have a graph similar to the one presented 
in Fig. [T3J as it is possible that E q2 weakly dominates strategy E qi and vice-versa. 
See [33J for details. 

The concept explained above clearly extend the concept of ESS for the PDE case. 
As the PDE case works as an approximation for large N of the discrete case, it is 
easy to see that we can extend the ESS definition also to the more realistic discrete 
case. 

Different to the most known ODE (see also the next section) case for the definition 
of ESS, here we cannot guarantee that the probability distribution will, in the long 
range (when adequately parametrized) accumulate in the ESS (when it is in the 
interior of the interval [0, 1]), but we can see that an individual that plays strategy 
I and II with frequencies given by the game's ESS is optimized to win any contest 
(with the same parameters). 

If the strategists involved in the game play with frequencies different from the ESS 
(for example, the pure strategies) the ODE prediction is that a stable mixture will 
evolve. This is impossible in the discrete case (as, in the long range, all individuals 
will descend of a single one in time t = 0, which will be of one of the given types) 
and also in the PDE model (as shown by Theorem |2J) . 

More generally, we say 

Theorem 4. Let PN,At( x 't) ^ e the solution of the finite population dynamics (of 
population N, time step At = 1/N 2 ), with initial conditions given by p° N (x) = p°(x), 
x = 0, 1/iV, 2/iV, • • • ,1, for p° e L\([0, 1]). Assume also that (A - 1, B - 1, C - 
1,D — 1) = 1/N(a,b,c,d) + 0(1/N 2 ). Let p{t,x) be the solution of the continuous 
model with initial condition given by p°(x). If we write p™ for the i-th component of 
PN,At( x it) i' n the n-th iteration, we have, for any t* > 0, that 

lim p'S =p(t,x), re G [0,1], te[0,t*]. 

iV^oo 

Proof. First, we consider the matrix M obtained from M by deleting the first and last 
rows and columns. Then, we observe that the derivation of the thermodynamical limit 
shows that the discrete iteration given by M is consistent — in the approximation 
sense (3U] — with Equation without any boundary conditions, provided that we 
set A = 1 + a/N, and similarly for B, C and D. From the results of Appendix A, we 
know that the discrete iteration is stable, since <r(M) C (—1, 1). From Appendix B, we 
see that the continuous problem without boundary conditions is well posed in the D s 
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spaces defined there. In this case, we can then invoke the Lax-Ricthmyer equivalence 
theorem [30 to guarantee that the discrete model converges to the continuum one, 
in the limit At, Ax — > 0, with At = (Ax) 2 . More precisely, the iteration defined by 
M converges to q(t,x), the smooth part of p(t, x); cf. appendix FBI 

Now returning to the iteration defined by M. In order to finish the proof, we 
only need to show that P(t, 0) and P(t, 1) converges weakly to the appropriate Dirac 
masses. We shall do the computation for x — 0, the case x = 1 being similar. 

For x = the iteration defined by M reads 

P(t + At, 0) = P(t, 0) + ip ft,± 
Thus, letting t = and solving the recursion, we have that 

171 — 1 



1 / 1 \ 

P(mAt, 0) = P(0, 0) + - J2 P /v • 

3=1 ^ ' 



3= 

Since ei G converges weakly to So as — > oo — by considering test functions with 
support contained in (1/N, 1/N) — we need only to show that it has the correct mass 
at each time t. For this, notice that 

V (jAt, 1 J = V (j At, 0) + ^d x V (jAt, 0) . 

Since p(t, x) = NV(t, x), we find that, in a weak sense, 

lim p^ 2 -> / g(s,0)ds + P(0,0). 



□ 



4 The diffusionless case and the replicator dynam- 
ics 

We shall see in this Section that the ODE Replicator dynamics is equivalent to the 
diffusionless version of Equation (|25p. This will have important consequences that 
we shall discuss later on. Notice also, cf. Remark El that this is the correct limiting 
equation, if the payoffs decay slowly to one as iV — > oo. 
Thus, we consider 



d t p = —d x — x)(x(a — c) + (1 — x) {b — d))p) 
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cf Q hi P 


nn nip 


a > (3 > 


X* < 


X = 1 


X = 


a > > f3 


x* e (o,i) 


X = and 1 


X = X* 


> a > (3 


X* > 1 


X = 


X = 1 


> (3 > a 


X* < 


X = 


X = 1 


(3 > > a 


x* e (o,i) 


X = X* 


X = and 1 


(3 > a > 


X* > 1 


X = 1 


X = 



Table 2: Stable and unstable equilibria in the range [0, 1] for the non-degenerated 
(a 7^ (3 7^ 7^ a) replicator dynamics 




Figure 14: Flux of the replicator equation for pure strategy dominated games (above), 
mixed strategy dominated (below, left) and bistable games (below, right). Here, 
X* = (3/(j3 - a). Compare with Figure ITHl 

A weak solution of this equation is given by p(t, x) = 5x(t), if X(i) solves 

X = X(l - X)(X(a - c) + (1 - X)(b - d)) , (28) 

which is the simplest replicator equation [IS] for the two-person game with pay-off 
matrix given by 

(: 5) • < 29 > 

The stationary points of Equation (j2Ej) are given by 0, 1 and X* := (3/ {[3 — a). 
The most interesting scenario occurs when a < (3 and X* e (0, 1) (i.e, (3 > > a): in 
this case the only stable equilibrium is the non-trivial X = X*. For the full analyze, 
see Table El Compare also with the description of dominance in the previous section. 

Our definition of dominance seems more general than many definitions that appear 
in the literature E3 El EH] ■ Furthermore, the use of the thermodynamical limit 
in the analysis make it much more simple to work. In particular, consider a game 
between E qi - and .^-strategists and a given replicator dynamics such that any non- 
trivial initial conditional converges in t —>■ oo to one of the two trivial equilibria, say, 
X = 0. The replicator dynamics is given by 



X = X(l - X)(X( gi - q 2 f{a -f3) + ( gi - q 2 ){q 2 a + (1 - q 2 )(3)) . (30) 
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If, for any initial condition X(0) G (0,1), lim^ 00 X(t) = 0, then, X(t) < 0, VI 6 
(0, 1), Vt G R+, i.e., (X( gi - g 2 ) 2 (a - /?) + ( 9l - g 2 )(g 2 a + (1 - q 2 )f3)) < 0. This 
implies that F^ qi q2 ^{y) > 0, G (0, 1), where F( quq2 ) is defined by Equation (|27ji. In 
particular F^ qim -\ is increasing and then from Lemma 121 we have .E g2 >- If, on 
the other hand, lim^ooX(t) = 1, by a similar argument, we have that E qi >- E q2 . 

These picture is completed after looking to Figures [T3] and |U and noting that the 
flow of the replicator dynamics always goes from the less-dominant strategy to the 
more dominant one, if we consider an equivalence (at the replicator dynamics level) 
between mixed populations of pure strategist and populations of mixed strategists. 

In reference (22] a thermodynamical limit of a frequency-dependent Moran process 
was also designed, but the pay-off were not re-scaled when N — > oo and the Fokker- 
Planck equation obtained was claimed to be valid for large, but finite N, and not in 
the thermodynamical limit. 



5 The Frequency Independent Moran Process 

In order to consider frequency-independent fitness, we impose a pay-off matrix such 
that the gain of a player is independent of others player's strategies, that is, A = B 
and C = D. In particular, we impose C/A = D/B = r. The number r is know as the 
relative fitness. Most results here are simple corollaries of results from the previous 
section. We state them only for completeness. 

Corollary 3. The fixation probabilities F n of type A individuals for an initial condi- 
tion of n mutants in the frequency independent Moran process with relative fitness r 
are given by 



1 — r 

n 



l-N 



+ 



k r n 



r 



l-n 



N 1 



,1-N ' 



F. = jf, r=l. 

Proof. When the relative fitness is constant, i.e. pAr(n) = 1/r, ([13)1 becomes 

(p -!)(*-!)' 



(31) 
(32) 



fc=i H v 
N 1 



k=l 



N - 1 
N - 1 



n -l 



(33) 



We sum the series and prove the corollary. If r = 1, it is straightforward to see 
that F n = n/N. 

□ 
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Remark 5. In the case of birth/death process we have instead: 

1 „ „ 1 — r k 1 — r k 

k=l r fc=l 



^ = G iE^r = G iE ^ = G - 1 _ r - 1 _ r ^ 



where we used that 



1 — r 



A' ' 



Note that the coefficients obtained in Corollary |3] are different from the one ob- 
tained in [21] , which are the same as in Remark The difference is the result of 
differences between a death/birth and birth/death processes. Anyhow, the formulas 
are equivalent for large N. 

We also define 7 := a = (3 and then Equation (J25j) is 

d t p = d\ (x(l — x)p) — rydx (x(l — x)p) . (34) 
As a simple consequence of Theorem |2] for Equation (J34)) . we have 

Corollary 4. Let p be a solution of Equation \3J$ with initial conditions p° G 
L\ fl L°°([0, 1]). Then, in a weak sense, p°° := lirm ; _ +0O p(-, t) = 7r [p°]5o + 7Ti[p°]5i. 
Furthermore, we have 7r [p°] = 1 — 717 [p ] and 

r 01 = I' Io^ X P °{x)dx 

1[P 1 1-e-T 
If we start with p° = 5 Xo , then 

1 — e~ lx ° 

/ 7 (x ) := tti^o] = - _ e _ 7 , (35) 

and lim 7 ^ ^i^xo] — ^o- This is true because the neutral case corresponds to 7 = 0. 
Note that / 7 (0) = 0, / 7 (1) = 1, V7 and that / 7 (x ) > x if and only if 7 > 0. So, in 
the language of previous sections, A y B •<=>- 7 > 0. 

It is important to compare the probability of fixation in the continuous limit, 
Equation (|3*3|l. and the result obtained for finite population, Equation (f3*T|). To un- 
derstand the idea we should consider that, in the finite case, we have initially a fixed 
proportion k G (0, 1) of mutants, such that the probability of fixation is given by 

„reiV j.kN „/tjV— 1 j j.kN 

I _ r Ar-i — K x — r^ -1 ~ 1 — r N ' 



when N is large and r close to 1. To be more precise, if r(N) = 1 + j/N, 
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In order to compare that formula with (|35|) for large N, we need only to impose 
k = Xq (the initial fraction of mutants) and then e~ 7 ~ r~ N , i.e., 7 ~ N(r — 1), for 
r - 1 < 1 (valid for large N), in agreement with 7 = liniTv^oo N(r — 1) (compare 

with A22D ). 

We cannot avoid the comparison of our result with the classical results by Kimura (20] . 
Following this reference, let u(t,y) be the probability that a mutant allele, initially 
with frequency y and relative fitness s be fixed after a time t in a randomly mating 
diploid population of size N . Then 

d * u = V -^W^ d l u + °VQ- ~ y)d y u . (36) 

4iVo 

This equation and Equation (J3%|) are associated backward/forward Kolmogorov 
equations with suitable rescalings [Tlj . Then, for example, Equation (|H5jl is the same 
found in [20 , where 7 = ANs for s = r(N)—l is the selective advantage. Furthermore, 
the fact that fo(x ) = x reproduces the idea that a neutral mutant (7 = 0) is fixed 
with probability equal to its initial frequency. 

Following, again, reference [Hj, if u(t,y) solves Equation (JHSJ), then u(t,x) = 
us{x)p(x, t) where 

I _ e -N sx 
U ^ = 1 _ e -47V 0S 

is the stationary solution of Equation and p(x, t) solves (with appropriate 
rescalings and normalizations). This shows the equivalence of this deduction and 
Kimura's one. 

6 The drift limit 

The "drift limit" means that the time-step is re-scaled according to At = 1/N. In this 
case, we do not need to consider the weak selection limit, i.e., pay-offs (and fitness) 
are considered time-step independent. This problem is mathematically well posed, 
but, as explained below, it seems not to be an interesting limit from the modeling 
point of view. We state it only for completeness. 

First, we see what happens for the drift limit of the frequency dependent Moran 
process, i.e., Equation (|2~^) . 

Theorem 5. Let p be the solution of Equation \2J$ with initial conditions given by 
p° E L\ R L°°([0, 1]). Then, p°° = tt 5o + 7T*8 X * + 7fi5i ; where 7Tq + 7f* + 7fi = 1 and 
x* = -(B - D)/(A - B - C + D). Furthermore, if A - C < 0, then vr = 0; if 
B - D > then irx = 0; and if {AD - BC)/((A - C)(B — D)) < then vr* = 0. // 

X* £ [0,1], 7f* = 0. 
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Proof. We multiply Equation (|24jl by 

= (1 - a;) A /^- c )x- D /( B - 1) )(a;(A -5-C + £>) + £?- D )(^-bc)/((a-c)( B - C) ) 
and integrate from to 1. Then 

S 'X ^ = "I - B _ c + J) + ,( B + c - 25) i)<fa 

"1 

i/j(x)p(x, t)dx . 

From Gronwall's inequality, we find that is supported at the zeros of ip(x). □ 

Suppose that we have a game where the strategy I dominates (e.g., the Prisoner's 
dilemma, where strategy I means "defect"), i.e., A > C and B > D. If AD — BC > 0, 
7f* = 0, and if AD — BC < 0, then x* > 1, and this implies 7f* = 0. Eventually, the 
full population will play strategy I. 

For A < C and B < D, the full population will play strategy II. 

For the Hawk-and-Dove game we have A — C < and B — D > 0. This implies 
that (AD - BC)/((A - C)(B -£>))> and then p°° = 6 X *, where x* G (0, 1). 

Finally, for coordination games, A — C > and B — D < 0, then (AD — BC) / ((A — 
C)(B — D)) < and p°° = 7fo^o + tTi^i- To obtain the values 7fj, i = 0, 1, note that 
x* G (0, 1) and 



d t / pdx = , d t pdx = 



This implies that 



— / p°°dx = I p°dx 



o 

7fi = / p°°dx = / p°dx . 

J X* J X* 

In a pictorial way, all the mass to the right of x* will move toward the point x — 1, 
while the mass on the left will move toward 0. If the initial condition is of delta-type, 
i.e., p° = 5 x o then the final condition is fully determined, p°° = 5q {p°° = Si) if xq < x* 
(xo > x*, respectively). 

Now, we consider the frequency independent case, i.e., we impose A = B = 1 and 
C = D = r at Equation (1241) . 

Corollary 5. Let p be the solution of 

x(l — x) 
x(r — 1) + 1 J 

with p° G L\ fl L°°([0, 1]). Then p°° = 5 1 for r > 1 and p°° = <5 /or r < 0. 



= ~(r - l)d x 



(37) 
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Proof. Note that ?p(x) = (1 — x) 1 ^ 1 r ^x r \ Then, its zeros are at most and 1. 
This implies 7f# = 0. The values of ttq and 7fi follow trivially. □ 

As a conclusion of this corollary, we note that the time-step of order 1/N implies 
in no diffusion, i.e., no genetic drift. So, the result of Equation (|H7j) is deterministic, 
in the sense that an arbitrarily small fraction of advantageous mutant will eventu- 
ally take over the entire population, while disadvantageous mutants will certainly 
be extinct (if the population is initially mixed). In Equation (|34|) nothing similar 
happens. 



7 Final Remarks 

The procedure used here can be applied to different evolution process. For example, 
consider the imitation dynamics given by the following rules: from a population with 
size A" and two possible types, we choose two individuals I\ and I 2 . If they are of the 
same type, nothing changes. If I\ is of type A and I2 of type B, I\ changes its type 
with probability ^(4>b — <Pa) and the same if we swap I\ and I2, where <pA and 4>b 
are the fitness for the types A and B respectively and ^ : R — ► [0, 1] is a continuously 
differentiable non decreasing function. Then, the transition coefficients are given by 

_/\T — Tl TL 

c+(n,N) = - N _^ !{(t> A -(t>B) , 

/ »rt 71 N — n _ . , ,. 
c.(n,N) = -__-^(0 B -0 A ) , 

c (n,N) = 1 - c + (n, N) — c_(n, N) . 

We consider the functions of x = n/N as defined in (|16|) - (jl8j) and with assump- 
tions fl2U)-(|22l we get 

lim N (cl ] + cj 1} + c a)N ) = 6x 2 *'(0)(a -b-c + d) 

+«'(0)(-a + 2b + c - 2d) - 2tf'(0)(6 - d) , 
Xim N [cf - c^A = -2x 3 ^'(0)(a -b-c + d) 

-2xV(0)(-a + 2b + c-2d) + 2x^'(0)(b - d) , 
lim (cf + cf + c m ) = -4^(0) , 

lim (c { l ] - c { lA = 2tf(0)(2x-l), 

lim (cf + c {0) ) = 2x{l - ar)tf (0) . 

AT— >oo V / 
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Gathering everything in Equation (|15p we find as the drift-diffusion limit of this 
process 

d t p = ^(0)d 2 x - x)p) - 29f'(0)d a {x(l - x)(xa + (1 - x)(3)p) , (38) 

with a = a — c and (3 = b — d. From the assumptions, ^(0), ^'(0) > 0. Relation 
of dominance for E gi - and E q2 -strategists are exactly the same as before, as can be 
easily computed from the fact that the conservation laws associated to Equation (|58j) 
are ip(x) = 1 and 

4>(?) = J exp ((-\(.Qi - Q2) 2 (a -(3)+ y{q t - q 2 )(q 2 a + (1 - q 2 )f3)^ dy . 

The coefficients can be adjusted from the basic discrete process. In particular, 
we can choose ^ such that Equation is drift-dominated (if ^(0) <C ^'(0)) or 
diffusion-dominated (if ^'(0) <C ^(0)). In a forthcoming paper, we will completely 
study this equation and this two different regimes. In particular, we can define a 
family of functions \l/ e , such that lim e ^ x I / e (0) = 0, but lim e _» ^e(0) > an( i use 
singular-perturbation theory to understand the diffusionless limit of the replicator- 
diffusion equation ([38j). We can expect a behavior similar to the one found in Sec- 
tion This means that, for certain imitation dynamics and for intermediate times, 
the evolution of the system, or more precisely, the "peak" of the density distribution, 
can be modeled by Equation (|25|l. as we can see in Figures ITU1 and ITT1 
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A Proof of Proposition U 

The following properties of /jv and will be useful in the sequel: 

1. M0) = f N (N) = 0; 

2. g N (0, r) = 1 - l/N and g N (N, r) = r - l/N; 

Notice that the first column of M is e± and last one is e/v+i- Hence 1 G cr(M), and 
ei,ejv+i are associated eigenvectors. Also, since M is nonnegative tridiagonal, we 
must have cr(M) C R. 

Since M is column stochastic, and since all diagonal elements are nonzero, an 
application of Gersgorin theorem to shows that er(M) C (—1,1]. 
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We now show that the 1 is an eigenvalue of multiplicity two. First, observe that 
M has the following block structure 




where M is a (N — 1) x (N — 1) tridiagonal matrix, with nonzero elements in the 
super and subdiagonal. Hence, M is irreducible. 

Let rji denote the sum of elements of the i-th column of M. Then we have 

r)i — 1, i = 2, . . . , N — 2 and < rji, r/^^i < 1. 

Because of the irreducibility of M, the strict inequality for 771 (or r/jv-i) is sufficient 
to show that 1 g <r(M) (cf. HE|). 

This result on the spectrum of M, together with the block structure of M proves 
the claim. 

We write 

M = PAP-\ 



where 



(\ 



P 



* * * 



0\ 



* * * 
yo * * * \J 



and A 



'I 0...0 N 
J 
v 0...0 1, 



We also notice that P -1 has the same structure of P. 

From the localization results on eigenvalues of M, we know that a (J) C (—1,1] 
and hence 



In this case, we have that: 



lim A* 

k^rOO 



lim J k = 0. 

fc^oo 



/1 



0\ 



V° ••• V 



and the result follows. 
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B Proof of Theorem U 



First, if we let a = a — c and (3 = b — d in ([25)1 . we have 

d t p = dl[x(l — x)p) — d x [x(l — x)((3 + (a — (3)x)p) 

Further, let 

e^ x+{a - 0) ^)w(t,x) = x(l - x)p(t,x). 
Equation (|3~9")l then becomes 

d t w — x(l — x) IdPw — — + j {(3 + (a — (3)x) 2 w 



(39) 



(40) 



First, we observe that by writing w e = w + et allows us to prove a maximum 
principle for C 2 (M + x (0, 1)) solutions to (|4U|) in a standard way. In particular, since 
w > is in the parabolic boundary, it is nonnegative everywhere. 

Existence can be established by Fourier series theory. In what follows, all the 
Banach spaces in this section are weighted with respect to 



oj(x) 



x(l — x) 



(41) 



Consider the associated equation 



if; = \uj(x)ip, ^(0) = ip(l) = 0. (42) 

The relevant facts 



Since w(x) G Ll oc ((0, 1)), standard Liouville theory applies to 
are collected in 

Lemma 4. Equation defines a singular Sturm- Liouville problem satisfying the 
following: 

1. The extreme points are singular points of limit point, non- oscillatory type. The 
Friedrich's extension of the operator on the left hand side of fcfify is a self-adjoint 
operator in H 2 ([0, 1]) PI -£^°([0, 1]), that is bounded from below. 

2. The eigenvalues of J7^ ) are real, purely discrete, bounded from below, and ac- 
cumulate only at infinity. 

3. The associated eigenfunctions are an orthonormal basis o/L 2 ([0, 1]). 
4- If denotes the spectrum, we have 

lim = K ^ 0. 
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Proof. A straightforward Frobenius analysis near and 1, shows that only one of the 
linear independent solutions can square integrable with respect to u(x). Moreover, the 
Frobenius expansion are regular without complex exponents. Hence, the extremes are 
of limit point, non-oscillatory type. The other results are standard — see for instance 

m □ 

An important property of (J42j) is given by 
Lemma 5. The operator defined by fl^jj| ) is positive- definite. 

Proof. For a > (3, this is straightforward. Also, since (f4*2*|) does not have continuous 
spectrum, the eigenvalues are continuous functions of the parameters. Hence, it is 
sufficient to show that zero is not an eigenvalue of ()42|) when a < (3. 
Thus, letting A = 0, and £ := (3 - a > in (HI yields 



~ (3 - - i 

4 2 



<p = o, p(o) = p(i) = o. 

which can be further transformed by letting x = £~ 1 f3 + v / 2£~ 1//2 ?/ in 

^ - (y 2 - l)(p = 0, <p(A)=<p(A + B) = 0, (43) 



where 



A = -— £~ 1/2 /3 and B = ^ 1/2 - 
2 2 



The general solution to (|43|) is given by 

<p(y) = e~ y2/2 ^ci + c 2 J e s2 ds 

On applying the boundary conditions, we see that a nontrivial solution exists if, and 
only if, we have 

= / e s2 ds - / e s2 ds = / e s2 ds. 

JO JO J A 

The last equality and the positiveness of the integrand implies B = 0, and hence 

e = o. □ 

Proposition 2. The initial value problem defined by Equation \40\) and w(0, x) = 
w (x), with w G L 1 ([0,1]) is well posed andw(t,x) G C°°(M + x [0,1]). Furthermore, 
we must have 

lim w(t, x) = 0, x G [0, 1]. 

t^oo 
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Proof. Let ipj satisfy (H2J) with Xj. Given / 6 L 2 ([0,1]) we set 

/ = £ /(7>i 
i>o 

Also, as in [31], define for s G R 

P s = (^£^([0,1]) 

Now, let 



5>(j)| 2 A*<oo 



w(t, x) = £ ^ (j>~' Aj ^)- (44) 
i>o 

For t > 0, it is clear that if satisfies (153)1 . If w G r> s for s > 1/2, then we have 
a classical solution. In any case, however, notice that (|44j) implies that w(t,x) G 
C°°(R + x [0,1]), and that 

lim w(t, x) = 0, 16 [0,1]. 

t— >oo 

□ 

Furthermore we have 

Lemma 6. Assume that wq G -£^([0, 1]) and let I{t) = w 2 (t, x) dx. Then, we have 

I(t) < I(0)e- 2Xat . 
Proof. From the Fourier representation of w(t,x), we have that 

oo oo 
3=0 j=0 

□ 

The solution given by (jB)) , while well defined and quite regular, has a major 
drawback: it does not satisfy, in general, the required conservation laws, as it can be 
checked by starting with a positive initial condition, and hence with positive mass. 
But the decaying property of the (|44|) implies that the mass will go to zero as time 
goes to infinity. 

We shall give up as little regularity as possible, and look for a solution in the class 
C°°(R + x (0, 1)). Thus, we shall write 

p(t, x) = q(t, x) + p D (t, x), (45) 
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where q(t,x) satisfies (|3*§|) without boundary conditions, and p-o(t,x) is a distribution 
solution with support in (0, oo) x {0, 1}. In this case, we must have, for some pair of 
nonnegative integers M and M' that 

M M' 
p D (t,x) = J2^(t)S k +Y,h(t)5 k 1 , (46) 

fc=0 k=0 

where 5^. means the k-th derivative of the delta distribution at Xn. 

Before proceeding, we must indicate precisely what we mean by a weak solution 
in this case. 

Definition 2. A weak solution to Actfy) will be a distribution with support in [0, 1] that 
satisfies 

"OO pi /*oo /*1 

/ p(t, x)d t 4>(t, x)dxdt = / / p(t, x) [x(l — x)dl<p(t, x) + 
o Jo Jo Jo 

+x(l -x)(p+(a- P)x)d x (f>(t, x)] dxdt+ 

+ / p°(x)(j)(0, x)dx, 
Jo 

where 

4>(t,x)eC™([o,oo) x [o,i]). 

Remark 6. Notice that the test functions in definition^ are required to be of compact 
support in [0, 1] and not just in (0, 1) as usual. Similar definitions have been given in 
other contexts; see for instance Wty . 

This definition can be recasted in the framework of usual distribution theory, by 
introducing the compactly supported distribution 

M M' 



u(t,x) = J2a k (t)5% + ^2b k (t)5'{ + X[o,i](x)q(t, 



k=0 k=0 

where X[o.i] is the characteristic function of unit interval. In this case, the distribution 
can act in C°°(M) and its entirely determined by the behavior in the support; see for 
instance JT^j/. We shall abuse language and shall, henceforth, identify u(t,x) with 
p(t,x). 

We now can state the following important result: 

Lemma 7. 1. Given p°(x) 6 L 1 ([0,1]) ; there is a unique weak solution p(t,x) of 
Efflj) such that p{t, x) e C°°(R + x (0,1)) that satisfies 

d f 1 d f 1 

— / p(t, x)dx = and — / ip(x)p(t, x)dx — 0. 
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2. This unique solution can be written as 

p(t, x) = q(t, x) + a(t)5 (x) + b(t)5 1 (x), 

with 

a (t) — / <?(s, 0)ds and b(t) — / q(s, l)ds, 
Jo Jo 

where q(t,x) is given by jjJW - 

Proof. We begin by substituting with p D given by (pltjj) Definition to 

obtain 

oo ( M M' 

-(!) fe+1 / ^a fc (t)^0(t,O) + ^& fc (t)^0(t,l)+ dt 

^° lfc=0 fe=0 J 

OO /"I 

q(t, x)d t 4>(t, x)dt = 

max(fc,2) 




JO 



= / ^a fc (t)(-l) fc £ ( 3 )di[x(l-x)}\ x=0 dt 3+2 Ht,0)dt+ 
Jo k=o j=o W 

poo M' max(fc,2) , 

+ 1 E^x- 1 )* E (y^w i -^u=i^' +2 ^ l )+^ 



fc=o j=o 

"OO /"l 

52 




JO 



q(t, x)x(l — x)d x (j)(t, x)dxdt+ 

M max(fc,3) 



k=0 j=0 ^ ' 

„oo A/' max(fc,3) , .* 

+ 1 E^x- 1 )* E ^J^[*(i-*)G9+^)]u^- i+ V(*,i)dt+ 



fc=o j=o 

OO fl 

q(t, x)x(l — x) (P + rjx)d x cj)(t, x)dxdt+ 




o jo 

i 

0/ 



+ / p (x)(j)(0, x)dx. 
Jo 

In the calculation above, we used that 

d x [x(l-x)]=0, i>2 and [a;(l - x){(5 + r}x)) = 0, i > 3, 
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Using that q is smooth, integrating by parts, and using ()39l yields the following: 

oo ( M M> \ 

(_l)*+i / J^ afc (t)^0(t,O) + ^6 fc (t)^V(t,l)+ \dt = 



DC 







fc=0 k=0 

[g(t,l)0(t,l) + g(t,O)0(t,O)]dt+ 

M max(fc,2) 



J>(t)(-l) fc £ M W [x(l -x)] U = o^ +2 0(t,O)dt+ 

- fc=0 j=0 \ ' 

,>oo A/' max(fc,2) . - 

+ / $>(t)(-l) fc £ (n^[x(l-x)]U^+V(*,l) + d« 



fc=0 j=0 
oo M max(fe,3) 



+ /_ $>,(i)(-l) fe £ (^)^[o;(l-a ; )(/5 + r ? a;)]U = oa^ + V(t,0)dt+ 



fc=0 j=0 
-oo M ' max(fc,3) , 

+ / $>(t)(-i) fc £ (; M^[x(i-x)(^ + ^)]U =1 ^ + V(*,i)dt. 

fc=o ,-=o VV 

First, we look at x = 0. Since the above must hold for any test function we must 
have, for k = 0,1, that 

OO POO 

a {t)d t (f)(t,0)dt= / q(t,0)(p(t,0)dt 
n Jo 



a l {t)d t d x <f){t,0)dt= / ^a i (t)(-l) i ^[x(l-x)(/3 + r / a;)]| :c= o9 :c 0(t,O)dt 



'OO 

I 

JO 

For 2 < k < M, we have 

POO 

(-l) fc+1 / a k (t)d t %<f>{t,0)dt 
Jo 

k 



3 



i=0 



^ Z=fc-2 V * / 

min(fc+2,M) _ _ 

+ jf E ««(*)(-!)' j J J^ fe - 1 )[a;(l-a;)(^ + ^)]U = o^(t,0)dt. 
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For k = M + 1, M + 2, we find: 
= 



= r Eai(0(-i)'f z ( ! 2) Vr (fc - 2) wi-^)]u=o^(t ! o)dt+ 

^0 i=fc-l V ' / 



For = M + 2, the relation above is identically zero but, for k — M + 1, we have 
that 



/•oo 

= (-l) M M/ a M (t)<9f +1 0M)dt 



Hence om(*) = 0. 

Considering k — M, yields 



I 



poo 

1) M+1 / a M (t)dtd!f<j>(t,0)dt = 
Jo 

((— 1) M_1 (M - l)a M -i(t) + (-1) M M(M - l)a M (0 + (-l) M M/3a M (*)) «9f 0(t, 0)d*. 

o 

Since a M (t) = 0, we have that a M -i(t) =0 as well. For k — M — 1, we have that 

/•oo 

(-1) M / a M -i(i)d t 6f "VMM* = 

f° ((-1) M ~ 2 (M - 2)a M _ 2 (t) - (— 1) M_1 (M - 1)(M - 2)a M _ 1 (f)) d? <j>(t, 0)dt 
Jo 

pOO 

+ / ((— 1) M_1 (M - l)/3a M -i(0 + (-1) M M(M - 1)(77 - /3)a M (*)) 9f0(t, 0)dt. 



Again,we have auif) = au-iij) = 0; thus a m -2(f) = 0. 

For 1 < k < M — 2, we have a linear relation involving Oj(t), i = k, . . . , k + 3 
(when = 1, we have i — 1, ... ,3). If three of them are zero, then the remaining one 
is also zero. Thus, starting with k = M — 2 and proceeding inductively, we find that 
a,k(t) = for k — 1, . . . , M. Therefore, only ao(t) can be nonzero. 

An analogous argument shows also that only bo(t) can be nonzero as well. We 
now drop the subscripts and determine their values. 

Integrating by parts, the corresponding relation for a(t), we obtain 



poo poo pt 

/ a{t)d t <f>(t,0) = - / / g(s,O)ds0(t,O) 
Jo Jo Jo 
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Hence 



a {t) = / q(s, 0)ds + do- 
Jo 



A similar calculation shows that 



b(t) = f q{s,l)ds + b , 
Jo 



It remains only to show that the conservation laws are satisfied. Substituting the 
found solution on them, we find 

a'(t) + b'(t) - q(t, 1) - q(t, 0) = and a'(t) - q(t, 1) = 

respectively, which are obviously satisfied. □ 
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